Critical interfaces of the Ashkin- Teller model at the parafermionic point. 
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. ABSTRACT 

flj ' We present an extensive study of interfaces defined in the Z^ spin lattice representation of 

Cj I the Ashkin- TeUer (AT) model. In particular, we numerically compute the fractal dimensions of 

I ' boundary and bulk interfaces at the Fateev-Zamolodchikov point. This point is a special point 

t^ I on the self-dual critical line of the AT model and it is described in the continuum limit by the Z4 

parafermionic theory. Extending on previous analytical and numerical studies |10| 112] , we point 
out the existence of three different values of fractal dimensions which characterize different kind 

C^ ■ of interfaces. We argue that this result may be related to the classification of primary operators 

of the parafermionic algebra. The scenario emerging from the studies presented here is expected 
to unveil general aspects of geometrical objects of critical AT model, and thus of c = 1 critical 

^ I theories in general. 
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1 Introduction and motivation 

The study of the scahng Umit of interfaces in systems at critically has been show to be a very fruitful 
field of investigation which provided deep insights in the comprehension of critical phenomena [T]. The 
richness of conformal symmetry in two dimension (2D) make the two-dimensional critical systems an ideal 
framework to study these issues. In particular, there is a variety of 2D critical models for which exact 
methods of conformal field theory (CFT), combined with an available Coulomb-gas representation [2], allow 
the exact computation of all geometrical exponents characterizing the fractal shape of critical interfaces |3] . 
Among these models, we mention for instance the critical percolation, the self-avoiding walks, the loop erased 
random walks, or again spin lattice models such as the Potts models [T]. These studies have benefited from 
a great amount of numerical work [3] supporting the proposed theoretical scenario. In general, the critical 
models whose geometrical properties are well understood, even if there are often no rigorous proofs, can be 
associated to the critical phases of a one-parameter family of statistical models, the 0{n) loop models, the 
parameter n representing the loop fugacity. A remarkable recent development came with the introduction 
of the Schramm-Loewner evolution (SLE) which constitutes a family of conformally invariant stochastic 
growth process characterized by a parameter k [S] . The SLE approach offers a conceptually new description 
of certain boundary interfaces defined in the 0{n) models. 

The critical points of 2D systems can be classified according to different CFT families. Each family 
is characterized by a given set of (infinite) symmetries and the associated chiral current algebra. The 
representation theory of these algebras is at the basis of CFT constructions. The critical phases of the 
0{n) models are described by the most simple of such CFT family, i.e. the one associated to the conformal 
symmetry alone and thus to the corresponding Virasoro algebra [6]. The connection between this family of 
CFTs and SLE approach has been fully understood [7]. However, there are other families of CFTs, the so 
called extended CFTs, which, besides the conformal symmetry, enjoy additional infinite symmetries. These 
theories describe universality classes which are different from the ones of the 0{n) models. A variety of 
statistical lattice models described by extended CFT have been introduced and studied since long time. 
Typically, these models are characterized by symmetries of some internal degrees of freedom, such as for 
instance the SU{2) spin rotational symmetry in the quantum 1 + 1 spin chains [8; . The role of the additional 
symmetries, in particular for the so called rational CFT (RCFT), is well understood for many aspects, as, 
for instance, the operator algebra of the primary operators, or the classification of the conformal boundary 
conditions. Nevertheless, despite all the recent activity and progress, the geometrical properties of critical 
interfaces defined in such extended CFT are in general not understood. Moreover, an SLE approach to 
describe extended CFT is not known. For this respect, the Z^ spin models offer an ideal laboratory to 
study these issues. These are lattice model of spins which take N values and interact via a nearest-neighbor 
potential which is invariant under a Zn cyclic permutation of the N spin states. The Zn spin models admit 
critical points, the so called Fateev-Zamolochikov (FZ) points, described by the parafermionic CFT which 
are extended RCFTs with Zn symmetry [5^. The cases with N = 2 and iV = 3 correspond respectively to 
the Ising and the three-states Potts model, which in turn are related to the critical phases of the 0{V2) 
and 0(-\/3) loop models. This is also manifest in the fact that the associated Z2 and Z3 parafermionic 
theories coincide with the c = 1/2 and c = 4/5 Virasoro minimal model M3 and M5 6,. The role of the Zjv 
symmetry becomes instead crucial for iV > 4 where the Virasoro algebra is not rich enough to describe the 
corresponding parafermionic theory. 

In the Zi\r spin lattice models, the boundary and bulk interfaces can be naturally defined and the role 
of the Z]\j internal degree of freedom is quite explicit in their definition. Moreover the Z^ spin models 
are simple models to be studied numerically. For A^ > 4, a description of the spin interfaces in terms of a 
low-energy effective field theory is not known, as it is the case for the interfaces defined in the 0{n) models. 
For this reason, the geometric description of the Zff spin models are for many aspects unknown. The study 
of the Zjsi spin models at the FZ point is thus expected to provide general deep insights on the geometrical 
description of extended CFTs. 

Numerical measurements of the fractal dimensions associated to the spin interfaces for Z4 and Z5 spin 
models at the FZ point were presented in [101 El- In [IH] we considered these Z4 and Z5 spin models on 
a bounded domain and we investigated the properties of a boundary interface related to certain boundary 
conditions. The numerical results were in the agreement with the theoretical predictions in |12) where the 
fractal dimension of this interface on the basis of the hypothesis of an SLE with an additional stochastic 



motion in the internal group of symmetry. This approach has been inspired by previous work on the connec- 
tion between SLE and CFT with superconformal symmetries |13) or with additional Lie-group symmetries, 

[mn]. 

In order to further investigate the geometrical properties of Zn spin models and the possible consistency 
with some proposed theoretical scenario, we studied systematically in [11] the bulk geometrical properties 
of spin and random cluster interfaces for the Z4 and Z^ models. These results clearly marked a difference in 
the behavior of these non local objects compared to the Ising or the three-states Potts model. Among these 
results there was the observation that the fractal dimension of certain spin bulk interfaces were different 
from the ones corresponding for the boundary interfaces in jlO| . 

The Z4 spin model is particularly interesting as it coincides with the Ashkin- Teller (AT) model. In [15] 
the fractal dimension of a boundary interface was investigated along the critical line of the AT. At the FZ 
point of the Z^ spin model, the value of the fractal dimension was found consistent with the value of the 
interface studied in [T2] and [TU]. 

In this paper we consider in great detail the spin cluster interfaces of the Z^ spin model at the FZ point by 
studying systematically different bulk and boundary interfaces. In order to interpret the numerical results, in 
particular in the light of their universal character, we discuss the classification of the Z4 conformal boundary 
conditions in terms of boundary spin configurations. 

The main result of this paper is the computation of three different values of fractal dimensions which may 
encode universal geometrical properties of the Z4 FZ point. Even if this results lacks of a clear theoretical 
explanation, we point out, on the basis of the properties of the Z4 CFT, a possible scenario for the geometrical 
properties of this extended CFT. It is important to stress that the Z4 CFT coincides with a free Gaussian 
field compactified on an orbifold. The study of the AT interfaces complement thus the results known for the 
free Gaussian field compactified on a circle [TU [T71 [TB] and may suggest an emergent general behavior for 
critical interfaces in c = 1 critical theories. 

2 The model 

We first define the Z4 spin model. On each site of a square lattice there is a spin which can take 4 values. 
Si = 1,- ■ ■ ,4. The Hamiltonian defining the model under consideration can be written as 

Hz,=-J2 ^'^■5-s, + BSs„s,±i + C(5s,,s,±2 , (1) 

where A,B,C are real nonnegative coefficients, Ssi,s- = 1 if S'i = Sj mod 4 and otherwise. Besides 
a constant irrelevant term, the interaction is described by two independent real parameters. The cyclic 
Z4 symmetry of the model (|T|) is completely manifest. One has to remark that the above Hamiltonian is 
invariant under a bigger symmetry than Z4, namely the dihedral group D4 symmetry acting on the spin 
degree of freedom. The Boltzmann weight corresponding to (|T|) reads: 



eM-HzJ = n 
where, by setting the normalization xo = 1 



TriSi — Sj) 
xq + 2x\ cos — - — h X2 cos7r(5i — Sj) 



(2) 



_ exp(A) - exp(C) _ exp(^) - 2exp(B) + exp(C) 
^\ -^ ^2 ~ ^ • (3) 

The ^"4 spin model can be considered as a generalisation of fourth-states Potts model, obtained by choosing 
x\ = X2 {B = C). The generalization to xi 7^ X2 consists in the possibility of having a non trivial weight 
between two non equal spins Si and Sj which depends on the difference between these spins, i.e\Si — Sj\. 

The Z4 spin model (jlj is a representation of the AT model [12]. An equivalent (and more standard) 
representation of the AT model is in terms of two coupled Ising models. In this Ising representation, on each 
site i of a square lattice one associates a pair of spins, denoted by Ui and r^, which takes two values, say 
up(-l-) and down(-). The Hamiltonian is defined by 

Hat = - ^ K{ai(Tj + TiTj) + K^^aiCJjTiTj . (4) 



In this representation the two parameters, K and ^^"4, correspond respectively to the usual Ising spin inter- 
action and to the 4-spins coupling between two Ising models. One can pass from the representation ^ to 
the ([U via the correspondence: 

Si = l ^ (Tj = +,Ti = + ; Si = 2^ ai= +,Ti = - 
S'i = 3 — > (Tj = -, Tj = - ; 5i = 4 — > (Tj = -, Tj = + ; 

from which one can derive the following relation between {K, K4) — > (xi, X2): 

exp(4i\:) = ■ ; exp(2X + 2X4) = ■ (5) 

1 - 2xi + X2 1 — 2:2 

The AT model on the square lattice is equivalent to the staggered six vertex model. It presents a rich phase 
diagram which has been very well studied [201 121] : in particular the phase diagram shows a critical line 
which is defined, in the Ising representation, by the self-dual condition sinh2iir = exp(— 2ivr4) and terminates 
at coth2iir2 = 2. By imposing the self-dual condition, the AT model can be solved by mapping it to the so 
called F model, a special case of a solvable six vertex model [21] . 

On the critical line one can identify three particular points : i) the fourth-states Potts model K = K4 
corresponding to xi = 0:2 = 1/3), ii) the case of two decoupled critical Ising models, corresponding to 
y/x2 = a;i = — 1 + \/2 (where K4^ = 0) and iii) the so called Fateev Zamolodchikov (FZ) point defined by 

FZ _ sin(j^) FZ _ ^ ^^"(t|) (a) 

sin(fl) sm({|) 

The FZ point has been shown to be completely integrable [9] . The continuum limit of the AT model on the 
critical line is described by a free Gaussian field with action S — j^ J dzdzd4)d(t> where the scalar field is 
compactified on a orbifold of radius r"'"'', i.e. cf) = (j) + 2Trr°^'' and (f> = —(f). The critical Potts, (Ising)^ and 
FZ point corresponds respectively to r°^^ = 2, ^/2 and r"'"^ — \/3. 

The Gaussian theory with r°^'^ = \f?> describes the FZ point and coincides with the Z4 parafermionic 
field theory [22]. In this paper we will mainly focus on the geometric critical properties of the FZ point. 

3 Classification of ^4 boundary states and their spin representa- 
tions 

In this section we discuss the problem of the classification of the boundary states for the Z^ parafermionic 
theory. In particular we are interested in the representations of such states in terms of spin configurations. 
The reason is that we want to study interfaces which are generated by imposing special boundary spin 
configurations. A classification of conformal boundary states in terms of spin configuration can then be used 
to identify interfaces whose measure is conformally invariant. 

The relation between certain spin configurations on the boundary and the conformal boundary states is 
known in the case of the Ising and three-states Potts model [321 IHl Hi] • This is in general not true for the Z4 
theory (and general Z^v theory, A'^ > 4), where only the spin representations of a small subset of boundary 
states has been explicitly discussed |26) . On the other hand, we mention that a complete characterization of 
parafermionic boundary states in terms of A-D-E lattice models degree of freedom has been accomplished 
in [57]. 

3.1 Rational CFTs: classification of conformal boundary conditions. 

In the following we briefly review the algebraic formulation of boundary states for RCFTs [23] HH] ■ A RCFT 
is characterized by a certain chiral algebra and the corresponding Hilbert space contains a finite number of 
the chiral algebra irrep. \j > which closes under operator product expansion. 

One of the main results of the boundary RCFT [251 US] is the bijection between boundary conformal 
states, which wc indicate as |j >, and the irrep. \j >. In particular, the states \i > can be expressed as: 

|z>=^^|j> (7) 



where the matrix Sij determines the modular transformation properties of the partition function defined on 
a cyUnder. For a RCFT the S matrix is in general known. The notation |0 > usually indicates the (trivial) 
identity representation. 

We define Zj-^ the partition function on a cylinder with boundary conditions \i > and \j > at the two 
ends. One has: 

%-E-i7A:.(g) (8) 

i 

where q is the modular parameter, the nl^ is the number of copies of the representation \i > occurring in 

the spectrum and Xiil) is the character of the representation |i >. A general boundary condition changing 
operator (b.c.c.) ipg-j produces a transition from the vacuum |0 > (related to the identity representation) 
and j boundary condition. The b.c.c. ^-r transforms in the representation j. In this case the cylinder 
partition function reduces to a single character as n-- = 1 and 

2o3 = Xj{q)- (9) 

The formulas ([8][9|) are extremely useful: given a certain lattice model, one can compute numerically the 
corresponding partitions ([HE]) for different boundary conditions on the cylinder. Then, by comparing these 
results with the boundary GET data, one can in principle associate the CFT boundary states to specific 
configurations of the degrees of freedom defining the lattice model. 

We specify the above equations for the Zjv theories. Reminiscent of the coset construction of the 
parafermionic theories, Z^r = SU{2)ff/U{l), we use the notation (see [20]) \l,m >, with / = 0,1/2,1,.. 
and m an integer to label the primaries of the Z^ theory. In particular the set of distinct principal represen- 
tations of the Zn theory is given by pairs |^, m > where I — 0, 1/2, 1.., N/2, m = —21, —21 + 2.., 2N —21 — 2 
with 2/ + m = mod 2. The representations \l,m > and \N/2 — l,k + m > have to be identified. The 
corresponding boundary states |^, m > are defined by: 



l,m ^l,m 



where 5; ^ is the modular transformation matrix 



,rm' ^ 2 ^,^„,^7^ ^.^ iT{2l + l) {2V + 1) 

'-" ^N{N + 2)' N- 



Sl^' = ,^-^nra'/N ^-^ ny.^^.jy.. -^ . , ^^^^ 



3.2 A known example: Three-states Potts model 

Before considering the Z4 spin model, we would like to review the results in [23l [24l [25] concerning the 
critical Z3 spin model, i.e. the critical three-states Potts model. In this case one can identify (almost) all 
the boundary states in terms of spin configurations with quite simply arguments. Still, the Z^ model is 
enough rich to show some general properties of the Ziq spin models which we will try to generalize to the 
case A^ = 4. 

The table of principal fields of the Z^ theory with the corresponding conformal dimension A is presented 
in Tab. [T] The fields VP*^ are the symmetry currents generating the Z3 chiral algebra. It is useful to stress 
that the fields listed in the above table are primaries of the Virasoro algebra but not of the parafermionic 
one. Indeed one can group the -Z3 Virasoro primary fields into two families, {/,'i'^,\l/~^} and {$^,£,$~^} 
which in turn correspond to the two representation modules of Z^ algebra. The fields in each module are 
thus connected one to the other by symmetry transformations, or, in other words, by acting with the modes 
of the '^^^ fields. This means that, taking into account (fTTj). the boundary states in ([T0|) transform under a 
Z3 rotation as: 

1070 >^ |o;2 >^ IM >^ |070 > 



|l/2,3>^ |1/2,1 >^ |l/2,-l >^ |l/2,3> . (12) 



Field 


A 


|Z,m > 


I 





|0,0> 


*i 


2/3 


|0,2> 


$-1 


2/3 


|0,4> 


$1 


1/15 


|1/2,1> 


e 


2/5 


|l/2,3> 


$-1 


1/15 


|l/2,-l> 



Table 1: Principal fields for the Z3 theory. 



This two distinct Z3 "orbifolds" can thus be directly related to the two already mentioned representation 
modules of the Z^ parafermionic algebra. 

Conformal boundary states in the spin representation 

In the three-states Potts model, where the spins can take the values 1,2 or 3, the following boundary 
conditions have been considered [2S1 [23] : 

• free: the spins can take the values 1, 2 or 3 with equal probability 

• fixed: the spins take the value 1 or 2 or 3. There are thus three fixed boundary conditions. 

• mixed: the spins can take with equal probability the value 1 or 2 (1 + 2), 2 or 3 (2 + 3) , 1 or 2 (1 + 2). 
Again there are three mixed boundary conditions. 



The transformation (|T2|) greatly constraints the possible boundary states identifications. Indeed, once 
one identify the boundary state |0,0 > with the fixed boundary condition (say when the spin are fixed to 
the value 1), the states |0, 2 > and |0, 4 > have to be associated to the other two fixed boundary conditions. 
Then, by observing that the free boundary conditions are invariant under a Z3 rotation, one is led to associate 
the states 1 1/2,1 >, |l/2,3 > and |l/2,— 1 > to the other mixed boundary conditions. Finally, taking into 
account that Z(i|i^2) = 2^(i|i+3) (for the Z3 symmetry), one can write: 



2^(l|l) = XI, 2(i|2)=X*i, -Z(i|3) = x^-i 
2^(l|l+2) = X*i = X*-i = 2(i|i+3), Z(i\2+3) = Xe 



(13) 



To make the connection with the minimal model M5 [6J with central charge 4/5 one has to take into 
account relations of the kind xi — X(i.i) +X(4,i) or Xs = X(i,2) +X(i.3) where X{r,s) is the Virasoro character 
of the operator (j)tr,s) in the minimal Kac table. The boundary operator ip^ transforming in the representation 
e generates then the boundary conditions (1|2 + 3) and thus the interface 5'Li?24/5 discussed in [3T] . 

One last remark: the six bulk operators which we have considered so far do not complete the set of 
primaries of the Z^ parafermionic theory. In general, the space of representation of a Zj^ theory includes 
also fields which are associated to the non-Abelian elements of the dihedral D^r group [3^. For the Z3 
theory, the two principal fields in this sector, which we indicate as Rq and i?i, have respectively dimensions 
1/8 and 1/40. The free boundary conditions have been associated with the representation Rq [35] 



-Z(l|/ree) — Xflo ' 

while the representation i?i has been associated to a new boundary condition: 

-2(l|„eto) — XBi- 

The physical interpretation of the new b.c. in terms of spin variables is actually unclear [25] . 



(14) 



(15) 



Field 


A 


\l,m > 


I 





|0,0> 


*i 


3/4 


|0,2> 


$-1 


3/4 


|0,4> 


^2 


1 


|0,6> 


$1/2 


1/16 


|l/2,-l> 




9/16 


|l/2,3> 


$-1/2 


1/16 


|1/2,1> 


£ 


9/16 


|l/2,5> 


$1 


1/12 


|1,2> 


e 


1/3 


|i,o> 



Table 2: Principal fields for the Z4 theory. 



The conformally invariant boundary states for the three-states Potts model are exhausted [33] by the 
eight boundary states (fT3 |) . (fT4l) and ([T5|) . Finally, one can check these identifications by verifying that the 
partition function associated to all possible combinations of such conditions are consistent with the fusion 
of the principal fields. For instance, one has 



2^(l+2|l+2) — XI 

The above relation is consistent with the fusion $^ x $^^ = 



Xe 



(16) 



I + e, associated to the spin boundary configu- 
ration (1 + 2|1|1 + 2) in the limit in which the two points where the boundary conditions change approach 
each other. 



3.3 Spin boundary states in the Z4 spin model 

The principal primary fields of the Z4 theory are shown in TabO In this table, the fields ^^^, ^^ generate, 
together with the identity /, the Z4 parafcrmionic algebra. Moreover, the fields ^^, together with e are the 
only neutral fields in the Z4 theory. In general the number of neutral operators is equal to the dimension of 
the phase space of the Z^ spin model. It is peculiar of the Zn theory with N even, A^ = 2n, that one of 
these neutral fields coincides with the current ^". 

Analogously to the case of the Z3 theory discussed above, one can group the above fields into three 
representation modules, {/, ^^, ^^^, ^^}, {<I>^/^,£ ,e ,(f>^^/^} and {<i>^,e}. Correspondingly, the boundary 
states will transform under a Z4 transformation as: 



1070 >-^ 1072 >^ IM >^ |o76 >^ 1070 > 

|1/2,1 >-> |l/2,3 >-> |l/2,5 >^ |l/2,-l >^ |1/2,1 > 

|T7o>^ IM >^ |T7o> . 



(17) 



Conformal boundary states in the spin representation 



In the Z4 model the spin can take the value 1,2,3 or 4. It is natural to consider the boundary spin 
configurations which generalize the Z3 spin model. Using the notations of the previous paragraph, we 



consider spin boundary configurations which are the natural extension of the one seen in the three-states 
Potts model: 



• 



• 



• 



free: 1 + 2 + 3 + 4. 
fixed:l, 2, 3 or 4. 

mixed: Three types of mixed boundary conditions: a) 1 + 2, 2 + 3, 3 + 4, 1 + 4, b)l + 3, 2 + 4 and c) 

1 + 2 + 3, 2 + 3 + 4,3 + 4 + 1, 1 + 4 + 2. 



Contrary to the case of the three-states Potts model, the Z4 transformations of the boundary states, 
consistent with (|17p and (ITU|) are not sufficient to identify all the boundary states. We have thus computed 
numerically the partitions ([9]). The numerical implementation is done with a transfer matrix. In order to 
compute ZiMB\ ~ xa^ it is convenient to work with an infinitely long strip of width L with boundary condi- 
tions A and B |23j . The corresponding dominant character xa is identified by computing the subdominant 

corrections in : 

ai (c/24-A) , , 

log(Z(^|B))^ao + -^+7r^^-^5 >-. (18) 

Here ag is associated to the bulk free energy while ai is a boundary term corresponding to the fixed boundary 
conditions A and B. We did the following identifications: 

-2(111) = XI 1 2^(l|2)=X*ii 2^(l|4) = X*-i 2^(l|3) = X*2 
2^(l|l+2) = X<I.i/2 = X$-i/2 = Z(i|i+4), -Z(i|2+3) = 2(i|3+4) = Xe' = Xe" ■ (19) 

The mixed boundary conditions of type a) have been discussed in [53] where the infra-red behaviour of a 
Zjv theory with free boundary conditions has been studied. In [5B1, the RG flow is characterized by a flow 
from the free boundary conditions to the four stable fixed points associated to the spins taking a fixed value 
on the boundary. These results generalize the boundary RG flow studied in [3S] . With a fine tuning of the 
boundary perturbation parameters [26], the free b.c. flows to the others four mixed boundary states of type 
a). _ _ 

The identifications of the spin boundary states associated to |1,0 > and |1, 1 > remain still ambiguous. 
Indeed our analysis could not discriminate between the possible identifications: 

2^(l|l+3) = X*i 2^(l|2+4) = Xe (20) 



or 



-2^(l|l+2+4) = X*i -2^(l|2+3+4) = Xe • (21) 

Below in the paper we will discuss the boundary conditions (1 + 2|3 + 4) which, as we will see, are related 
to the interface studied in [15 . We have verified numerically that -E(i+2|3-i-4) — Xe- This is consistent with 
the identifications (|T51) and with the operator fusion $^/^e — e + .., determined by the operator algebra 
[30) . Note that the identification of the b.c.c. operator transforming in the e representation as the one which 
generates the condition (1|2 + 3 + 4) has been proposed in [12]. In particular this would support the fact 
that the b.c.c operator transforming as e is related to such boundary conditions |15| . 

One last remark: analogously to the Z3 case, besides the fields shown in Tab. [21 related to the abelian 
Z/i sector of the theory, there are a set of fields R related to the Z2 reflections of the dihedral Z4 group. We 
have verified that Zijree = X-Roj where Rq is the (twist) field with dimension 3/24 considered in [12] . 

4 Boundary Interfaces. 

We will first present the results obtained by measuring interfaces connecting two points of a lattice in analogy 
with the chordal SLE interfaces defined in the critical 0{n) model. As already discussed in section 2, the 
Z4 spin model generalize the fourth-states Potts model. In the Potts model, one can in general consider 
either geometrical interfaces bounding spin clusters |31) , i. e. group of spins with the same value, or Fortuin 
Kastelyn (FK) clusters [Ml |3S] . In [TT] , we explained that only geometrical interfaces can be defined for the 
Zm spin models since the FK clusters do not percolate at the critical point. 



We generated these geometric interfaces by simulating finite square lattices of size L x L with certain 
boundary conditions {Ai + A2 ■ ■ ■ \Bi + B2 ■ ■ ■)■ This notation means that we set one half of the boundary 
spins to take the values Ai + A2 + .. and the other half the values Bi + B2 + ■.■ Moreover, we impose that 
the change from one condition to the other is on the middle of the two opposite borders of the square lattice. 
Then, for each spin configuration, there is an interface defined as the line on the dual lattice separating the 
Ai spins connected to one boundary from the Bi spins connected to the other boundary. 

The result for the interface associated to the condition (1|2 + 3 + 4) was already presented in [10] where 
we found a value of the fractal dimension di compatible with the value 1 + 10/24 predicted in [12l. Note 
that, as we mentioned in the previous section, the scenario proposed in [12j is based on the hypothesis that 
the b.c.c operator associated to (1|2 + 3 + 4) transforms as in the representation |1,0 > of dimension 1/3. 
Naturally, this would be consistent with the identification (pij) . In order to have a more general picture, we 
have systematically analyzed boundary interfaces associated to differed boundary conditions, mainly inspired 
by the classification discussed in the previous section. 

Before presenting the main results, we briefly explain how the fractal dimension of these interfaces has 
been obtained. In order to generate the configurations, we choose the Ising representation (|4]) for which it is 
possible to use a cluster algorithm [36] much more efficient than standard Monte Carlo [11]. For each type 
of boundary conditions, measurements were performed on systems up to the linear size 1280. For each sizes 
we first determined the autocorrelation time t{L) and we average over N{L) x r(L) with N{L) — 1 000 000 
for L < 160, iV(160) = 500 000, A^(320) = 250 000, iV(640) = 100 000 and A^(1280) = 50 000. A typical 
value for the autocorrelation time is r(L = 640) ~ 10000 for the boundary condition (1|2) and of the same 
order for other boundary conditions and the same linear size. The fractal dimensions are obtained in the 
following way. For each linear size L, we determine the average length of the interface 1{L) which should 
scale as 

1{L) ~ L'^f . (22) 

Then we obtain the effective dimension as 



L1+L2 



\og{l{L^)/l{L2)) 
\og{Li/L2) 



(23) 



where Li and L2 are two different linear sizes of the square lattice. 

In Fig.[T^) we present the effective fractal dimensions so computed for the four type of interfaces associated 
to the (l|2 + 3 + 4), (l + 3|2 + 4) and (l + 2|3 + 4)and (1|2) boundary conditions. The conditions (l + 2|3 + 4) 
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Figure 1: Effective exponents obtained for different interfaces. The three dashed lines correspond to 1 + 
10/24, 1 + 7/16, 3/2 as discussed in text. 



or (1 + 3|2 + 4) are particularly interesting as they are, besides the condition (1|234), the only other ways 
for defining a single interface for this model. As one can see from Fig. (1^), we measured a value for the 



fractal dimension of the interface (1 + 2|3 + 4) which is shghtly bigger than the one measured for the (1|234) 
interface. Even with the large amount of statistics that we have accumulated for this point, it is not possible 
to decide if the asymptotic value will converge to the same value as for the condition (1|2 + 3 + 4) or not. 
Note that the condition (1 + 2|3 + 4) was also considered in the recent work who obtained similar results, 
namely df = 1.4226(13) [15l|37]. 

On the contrary, for the condition (1 + 3|2 + 4) one clearly sees that the value is completely different and 
seems to converge to the value ^2 — 3/2. Few comments are in order. For this condition we do not have 
data for larger sizes. This is due to the fact that the cluster algorithm is much less efficient for studying 
this interface. To be more precise, the cluster algorithm is not able to take into account the (1 + 3|2 + 4) 
boundary conditions since the algorithm works for each copy of the Ising model separately [36] . Thus one 
can only change a spin S = l{a = +;r = +) in either S = 2(ct = +; r = — ) by updating the copy a or in 
S = 4(a = — ;r = +) by updating the copy r. The direct change S ~ 1 ^i' S = 3 can not be done within 
this algorithm. Since on one border the spins S — 1 or 3, then these spins will be frozen, the same being 
true for the other part of the border where we impose S" = 2 or 4. Thus for this type of boundary conditions, 
we consider a more complicated algorithm in which we alternate cluster update with standard Monte Carlo 
updates on the border. Another intriguing observation concerns the obtained value ^2 — 3/2. In fact, for 
this choice of boundary conditions, a fractal dimension compatible with ^2 — 3/2 is obtained all along the 
AT line, [38 . Moreover, in the next section, we show that we find bulk interfaces with fractal dimension 
d2- The fact that we obtain the value ^2 appear to us natural, even if we do not have any strong argument. 
Indeed it is well known that the value ^2 is the fractal dimension of certain interfaces related to the level 
lines of the Gaussian free field compactified on a circle [16l [ITl [18] . As we have seen in the section 2, the Z4 
spin model on the critical line is also described by a c = 1 CFT which is the free Gaussian field compactified 
on the Z2 orbifold, -\/2 < r°^^ < 2. It is thus natural to ask whether these interfaces (1 + 3|2 + 4) can be 
somehow related to an SLE process with k — A. 

We also report here other cases of boundary conditions. They are (1|2), (1|3), (l|2 + 3) and (l|2 + 4). The 
cases (1|2) and (1|3) generate two interfaces: an interface separates the spins connected to the boundary "1" 
on one side and spins 2,3,4 on the other side while the second interface separates spins 2,3,4 with spins 
connected to the boundary 2 or 3. By symmetry arguments, one can see that each of these two interfaces is 
equivalent [1|. In Fig.[T^), we show the effective fractal dimension for the case (1|2) and we average other the 
two interfaces. The result for the case (1|3) is identical except finite size corrections for the smaller sizes, 
thus we do not show the fractal dimension for this case. The value obtained for the largest sizes is df ~ 1.435 
which is clearly different from di. We observe that this value, that we will call ^3, is very close to the value 
1 + 7/16. In [39l Uni E], it was shown that the interfaces related to an holomorphic operator with spin s 
are described by SLEi^is) with k^^^ = 8/(1 + s). Here we just notice that the value 1 + 7/16 corresponds 
to 1 + k/8 with k(*) — 16/k^'*^ and s — 3/4 which is the dimension (spin) of the holomorphic current ^^^. 
As seen in section 3, the parafermionic current ^^^ is associated to the (1|2) boundary conditions. This 
value is also compatible with the value reported for bulk fractal dimension in [llj . We will come back on 
this point in the next section. The last two cases are a little bit more complicated. For (1|2 + 3) we have 
again two interfaces but now these two interfaces do not need to be equivalent. The first interface separates 
1 from 234 while the second interface separates 23 from 14. In Fig. [TJd), we show the fractal dimension for 
each interface. It is clear that they are not equal. For both interfaces, we observe strong finite size effects, 
but the difference between the two fractal dimensions remains near constant as we increase the linear size 
L. And it seems that the fractal dimension associated with interface bounding the domain connected to 
the boundary with Si = 1 converges towards d^ while the fractal dimension of the other interface converges 
towards di. For the last case (1|2 + 4) we observe again two fractal dimensions. The one associated to the 
interface bounding the domain connected to the boundary with 5,; = 1 converges towards d^ again, while 
the second fractal dimension converges towards d2. 

To summarize the numerical findings, we obtain fractal dimensions which can be grouped in three parts. 



^In [3r similar boundary conditions were considered for the three-states Potts model and called "fixed". In this study, the 
authors separate the interfaces in a "composite" part and a "split" part and claimed that the fractal dimension associated to 
the "split" part (~ 1.589) was much larger than the one expected for spin cluster boundaries and obtained for "fiuctuating" 
boundary conditions (1 + 10/24 ~ 1.41667). In fact the difference in these fractal dimensions is due to the existence of strong 
finite sizes corrections for the " split" part with " fixed" boundary conditions as we have checked for the three-states Potts model 
and for the Z4 spin model 38 . 



A first set of fractal dimensions converge to a value close to di ~ 1 + 10/24 and are associated to either 
(l|2 + 3 + 4) and one of the fractal dimension of (l|2 + 3). A second set of fractal dimensions converge towards 
d2 — 3/2 and are associated to (1 + 3|2 + 4) and one of the fractal dimension of (1|2 + 4). The last set of 
fractal dimensions converge to a value close to d^ ~ 1 + 7/16 and are associated to the fractal dimension of 
(1|2) or the second fractal dimension for (1|2 + 3) and (1|2 + 4). 

Finally the fractal dimension for (1 + 2|3 + 4) takes a value between di and ds, it is difficult to conclude 
definitely for this case. 

5 Bulk interfaces 

In [llj we determined the fractal dimension of interfaces around finite clusters in the bulk. This measurement 
was done by considering any type of clusters, i.e. a geometrical cluster of spins of any fixed value surrounded 
by spins with a different value. But since we observed in the previous section that the fractal dimension 
can depend on the type of interface and on particular in the number of allowed values on each side, we will 
check now if a similar property also occur for finite size clusters. The goal is to check if the fractal dimension 
can depend on the number of allowed values of the geometrical clusters and how. In our simulations, we 
considered a square lattice with periodic boundary conditions. For a given geometrical cluster we computed 
the average area A defined as the number of spins inside a contour of length I. These two quantities are 
related to the fractal dimension in the following way : 

A{1) = f (24) 

with 5 = 2/(i/. To obtain this last relation, one note that A{1) = R^ with R the radius of gyration which is 
related to the length of the contour by I = R'^f . 

Wc can then perform a direct measurement of the distribution of the clusters in function of their length 
recording for each length the average area. In order to compare the different type of clusters, we first use the 
global Z4 symmetry for each configuration. We rotate the spins such that the majority of them take the value 
" 1" . Next we compute all the type of spin clusters. The simplest case is the clusters "1" which corresponds to 
spins of value "1" surrounded by spins taking another value. We can then define the associated distribution 
Di (A, I) which counts the number of such clusters with area A and length of the surrounding interface I. One 
can define three one-spin distributions, Di,D2,D3 (by symmetry, the distribution D4 is equivalent to D2). 
Next we have to consider two-spins distributions corresponding to clusters of spins taking two fixed values 
surrounded by spins taking the remaining values. We have to consider D12, D13, D23, D2i and similarly for 
the three-spins, we have to consider Z5123, £'124, D234. For each of these distributions, one can then extract 
a fractal dimension by using eq. (l24p . Finally in order to get a precise measurement of S, we compute the 
integrated quantity corresponding to 

Xilma.) = J '""' Ail) C Cat • (25) 

The values of d/ — 1/6 is obtained from the ratio X{lmax)/X{lmax/'2) for increasing Imax for the different 
types of clusters. 

A first result is that the fractal dimension is the same for any type of cluster with equal sign even if we 
have introduce an explicit breaking of symmetry by imposing the majority rule. We will only consider the 
case of D2 in the following. The same result is also obtained for the clusters with spins restricted to three 
values for which we will only consider the case of £'234- For the clusters with spins restricted to two values, 
we observe two behaviors, one for D12 with the same result as for D23 and a second one for D13 with the 
same result as for D24. 

In Fig. [21 we present the obtained fractal dimensions versus Imax- This figure contains many similarities 
with Fig.[T^). Again we obtain a fractal dimension 1.5 associated with clusters containing spins of value 1, 3 
or 2, 4. This fractal dimension is in very good agreement from the one obtained for the interface (1 -I- 3|2 -f 4). 
The second fractal dimension is associated to clusters with one value ((2)) and seems to be in correspondence 
with the interface (1|2). It is also natural to associate the fractal dimension for clusters with three values 
((2 + 3 + 4)) with the one from the interface (1|2 -t- 3 + 4). The last one (1 + 2) or (3 + 4) is then associated 
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Figure 2: Effective exponents in the bulk obtained from eq. (p5|) for Z4. 

to the fractal dimension of the interface (1 + 2|3 + 4). The interesting observation is that, in the bulk, the 
fractal dimensions for the clusters with three values seem to coincide with the one for (1 + 2) in the large size 
limit. This would then confirm that the two fractal dimensions for the interfaces (l|2 + 3 + 4) and (l + 2|3 + 4) 
would coincide in the large size limit. 

An important observation is that the fractal dimension ^2 — 3/2 is associated to objects arising naturally 
in the high temperature expansion of the AT Model [5J [3TJ |31]. Indeed, by replacing in ^ (o'i,Ti) by 
{(Ti,ti = UiTi) then clusters of either t^ = +1 or t^ = — 1 are natural objects to consider. Indeed, as shown 
in [21 1211 142] , after performing an Ising high-temperature expansion, the partition function defined with the 
Hamiltonian Q represses as 

Zh^t^ Y^ (tanh2ir)'+^ (26) 

graphs 

with the sum running on two types of (non intersecting) graphs : polygons C on the lattice coming from the 
expansion of ai and polygons T) on the dual lattice coming from the low-temperature expansion of ti and 
the total numbers of bonds on each lattice is I and d respectively. The boundary of the clusters forming 2? 
corresponds to the bulk interface with the fractal dimension d2. 

6 Summary and conclusions 

In this paper we have considered the Z4 spin lattice model at the FZ critical point. The FZ point is described 
in the continuum limit by the parafermionic Z4 CFT which is a rational CFT with extended Z4 symmetry. 
Contrary to the case of the CFTs based on Virasoro algebra, the behavior of geometric objects for general 
extended RCFTs is far to be understood. In order to provide new insights into this problem, we have 
considered different types of boundary and bulk spin cluster interfaces which can be naturally defined in the 
Z4 spin lattice model. In particular, we have computed their fractal dimensions at the FZ point. The basic 
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idea behind the analysis presented here and in previous works [TU1[TT1[T2] is that, analogously to the critical 
0{n) models described by Virasoro CFT's, the geometrical properties of a Z4 spin lattice model at the FZ 
point should be related to the classification and properties of the primary operators of the corresponding 
extended algebra, in this case the parafermionic algebra. 

First we considered the lattice model on a bounded domain. We studied interfaces which origin and 
terminate at boundary points and which are generated by imposing certain boundary spin configurations. 
We examined all kind of boundary spin configurations whose definition is naturally suggested by the Z4 
spin symmetry of the model. Despite the great number of interfaces we considered, the numerical results 
indicate that there are only three values of fractal dimension, di,d2 and d^, characterizing the critical 
behavior of these boundary interfaces. We showed that this result can be somehow understood on the 
basis on the classification of conformal invariant boundary states. For this purpose, we have discussed the 
identification of conformally boundary states, as predicted by the boundary RCFT, in terms of boundary 
spin configurations. We provided evidences that the existence of this three different values can be related to 
the different representation modules of the parafermionic algebra. For instance, the value di ~ 1 + 10/24, 
already measured in |10[ lll| . is associated to the b.c.c. operator of dimension 1/3 and compatible with the 
value predicted in [12]. Other interfaces, like the one studied in [15], which seem to have the same fractal 
dimension di, are associated to conformally boundary conditions which in turn are connected by fusion to 
this operator. This was in general checked on the basis of our identification of conformal spin boundary 
states and fusion rules of the correspondent Z4 primaries. For the values ^2 and ^3 there are no theoretical 
arguments to derive them. However we stressed that i) the value d2 obtained is very close to the fractal 
dimension 3/2 of SLE interfaces defined as level lines of free Gaussian fields and related to a b.c.c. operator 
of dimension 1/4 [16l [ITl [H]. This is particularly interesting as the Z4 is also a c = 1 theory which can be 
described by a free Gaussian field compactified on a orbifold; ii)the value da ~ 1 + (1 + s)/4 where s = 3/4 
is the dimension of the -E4 holomorphic current. It is natural to ask whether there is some connection to the 
S'L£'g/(i_(_s) interfaces which are related to an holomorphic operator with spin s [MlHUlBTj . 

We have further investigated different interfaces which has been opportunely defined in the bulk. The 
results obtained for the bulk interfaces confirmed the scenario emerging from the study of boundary interfaces. 
Indeed, we could observe there are indeed three different values of fractal dimension which are compatible 
with the values di , d2 and da . Interestingly the interfaces having the fractal dimension d2 correspond to the 
spin cluster boundaries which appear in the high-temperature expansion of the AT model [H HU |35] . 

We conclude by emphasising that our results on critical interfaces at the parafermionic point of the AT 
are expected to unveil general aspects of geometrical objects of critical AT model, and thus of c = 1 critical 
theories [38]. Therefore, as a theoretical understanding is still lacking, we believe that the results presented 
here can represent good motivation for studying critical interfaces in free Gaussian field on an orbifold. 
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